Quantum anomalous Hall states in the p-orbital honeycomb optical lattices 
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We study the quantum anomalous Hall states in the p-orbital bands of the honeycomb optical 
lattices loaded with the single component fermions. Such an effect has not been realized in both 
condensed matter and cold atom systems yet. By applying the available experimental technique 
developed by Gemelke et al. to rotate each lattice site around its own center [lHsj], the band 
structures become topologically non-trivial. At a certain rotation angular velocity f2, a flat band 
structure appears with localized eigenstates carrying chiral current moments. With imposing the 
soft confining potential, the density profile exhibits a wedding-cake shaped distribution with insu- 
lating plateaus at commensurate fillings. Moreover, the inhomogeneous confining potential induces 
dissipationless circulation currents whose magnitudes and chiralities vary with the distance from 
the trap center. In the insulating regions the Hall conductances are quantized, and in the metallic 
regions the directions and magnitudes of chiral currents cannot be described by the usual local- 
density-approximation. The quantum anomalous Hall effects are robust at temperature scales small 
compared to band gaps, which increases the feasibility of experimental realizations. 

PACS numbers: 03.75.Ss, 05.50.+q, 73.43.-f, 73.43.Nq 



I. INTRODUCTION 

The anomalous Hall effect appears in ferromagnets in 
the absence of external magnetic fields which was discov- 
ered soon after the Hall effect. The mechanism of the 
anomalous Hall effect has been debated for a long time, 
including theories of the anomalous velocity from the in- 
terband matrix elements [3] , the screw scattering Q , and 
the side jump [f|. Recently, a new perspective has been 
developed from topological band propert ies, the Berry 
curvature of the Bloch wave eigenstates |7H12| , which has 
been very successful. The Berry curvatures serve as an 
effective magnetic field in the crystal momentum space, 
leading to an anomalous transversal velocity of electrons 
when electric fields are applied [l(|. The anomalous 
transversal velocity of electrons gives the intrinsic con- 
tribution to the observed anomalous Hall conductivity in 
ferromagnetic semiconductors. 

The integer quantum Hall effect (QHE) is the quan- 
tized version of the Hall effect, in which Hall conduc- 
tances are precisely quantized at integer values. This ef- 
fect arises in the two dimensional electron gases in mag- 
netic fields with integer fillings of Landau levels. The 
quantization of the Hall conductance is protected by the 
non-trivial band structure topology characterized by the 
Thouless-Kohmoto-Nightingale-den Nijs (TKNN) num- 
ber, or the Chern number (l3l. H3|. 

In order to achieve a non-zero Chern number pattern, 
time-reversal symmetry needs to be broken, but Landau 
levels are not necessary. Integer QHEs can appear as a 
result of the parity anomaly of the 2D Dirac fermions 
[IB - [l7l ]. Haldane constructed a tight binding model in 
the honeycomb lattice with Bloch wave band structures, 
and showed that it exhibits quantum Hall states with 
v = ±1 [T7I ]. This effect is termed "quantum anomalous 
Hall effect" (QAHE) because the net magnetic flux is 



zero in each unit cell and there are no Landau levels. 
The Haldane model has been taken as a prototype model 
for QAHEs. 

The Hall effect has been generalized into electron sys- 
tems with spin degrees of freedom as "spin Hall effect" , in 
which transverse spin currents instead of charge currents 
arc induced by electric fields fl8l - [23| . Different from the 
Hall effect, the spin Hall effect maintains time- reversal 
symmetry. Topological insulators are the quantum ver- 
sion of the spin Hall systems, which exist in both 2D and 
3D systems. Their band structures are characterized by 
the /^-topological index [241-32]. These states have ro- 
bust gapless helical edge modes with odd number of edge 
channels in 2D systems [26j, l33|, |34| , and odd number of 
surface Dirac cones in 3D systems [30rl32l ]. Topological 
insulators have been experimentally observed in 2D quan- 
tum wells through transport measurements [35I ] , and also 
in 3D systems of B^Sbi-^, Bi2Te3, Bi2Se3, and Sb2Te3 
through the angle-resolved photoemission spectroscopy 
[36H39I ] and the absence of backscattering from scanning 
tunnelling microscopy spectroscopy [40l - t42j |. 

Among all these Hall effects mentioned above, only 
the QAHE has not been experimentally observed yet. 
Several proposals have been suggested to realize this 
novel Hall effect in semiconductor systems with topo- 
logical band structures by breaking time-reversal sym- 
metry, such as ferromagnetic ordering 43 46]. Because 
no external magnetic fields are involved, QAHE states 
are expected to realize the dissipationless charge trans- 
port with much less stringent conditions than those of the 
quantum Hall effect. This is essential for future device 
applications. 

On the other hand, the development of cold atom 
physics has provided a new opportunity for the study 
of QHEs and QAHEs. Several methods to realize these 
effects have been proposed including globally rotating 
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traps and optical lattices, or introducing: effective gauge 
potential generated by laser beams |47H54j . In particular, 
the Haldane-like models were proposed in Refs [52l - l54l |. 
Furthermore, the realization of the quantum spin Hall 
systems has also been proposed [55(. All these proposals 
involve experimental techniques to be developed. 

In a previous paper [5(| . one of us has proposed to 
realize the QAHE in the p-orbital bands in the honey- 
comb optical lattices through orbital angular momentum 
polarizations. This can be achieved by rotating each op- 
tical site around its own center, but there is no overall 
lattice rotation. The net effect of this type of rotation 
is the "orbital Zeeman effect", which breaks the degen- 
eracy of the onsite p x ± ip y orbitals. This gives rise to 
non-trivial topological band structures, and provides a 
natural way to realize the Haldane model. Increasing the 
rotation angular velocity induces the topological phase 
transition by changing the band structure Chern num- 
bers. In the regime of large rotation angular velocities, 
the band structures reduce into two copies of Haldane's 
model for each of the p x ± ip y orbitals, respectively. 

The main advantage of this proposal is that all the ex- 
perimental techniques involved are available. The hon- 
eycomb optical lattice was constructed long time ago 
[57| . Recently, the superfluid-Mott insulator phase tran- 
sitions of bosons have been observed in the honeycomb 
lattice by Sengstock's group [58] . The rotation technique 
has been developed by Gemelke, Sarajlic, and Chu 
They have applied it to rotate the triangular lattice filled 
with bosons to study the fractional quantum Hall physics 
d, For the purpose of studying QAHE, we only need 
to apply this technique to the honeycomb lattice and load 
it with fermions. 

This proposal brings a natural connection between the 
QAHE and orbital physics in optical lattices. Orbital is a 
degree of freedom independent of charge and spin, which 
was originally investigated in solid state systems. It plays 
an important role in superconductivity, magnetism, and 
transport properties in transition metal oxides. The key 
features of orbital physics are orbital degeneracy and spa- 
tial anisotropy. Optical lattices bring new features to 
orbital physics which are not easily accessible in solid 
state orbital systems. First, optical lattices are rigid and 
free from Jahn- Teller distortions, thus orbital degener- 
acy is robust. Second, the metastable bosons pumped 
into high orbital bands exhibit novel superfluidity beyond 
Feynman's "no- node" theory [3^, l59l - [63j , which does not 
appear in 4 He and the previous study of cold bosons. Ex- 
citingly, this unconventional type of BECs have been ex- 
perimentally observed recently [6~H.l65l]. Third, p-orbitals 
have a stronger spatial anisotropy than that of d and /- 
orbitals, while correlation effects in p-orbital solid state 
systems (e.g. semiconductors) are not that strong. In 
contrast, interaction strength in optical lattices is tun- 
able. We can integrate strong correlation effects with 
strong spatial anisotropy more closely than ever in p- 
orbital optical lattice systems (66l - [69j . Recently, we also 
extend the research of orbital physics with cold atoms 



into unconventional Cooper pairings, which includes the 
/-wave Cooper pairing [70( in the honeycomb lattice, and 
the "frustrated Cooper pairing" in the triangular lattice 

a 

This paper is as an expanded version of the previous 
publication of Ref. 56] on QAHE in the p-orbital band 
in optical lattices. We will also present new results in- 
cluding the chiral flat band structures which occur at 
an intermediate rotating angular velocity. The effects 
of the confining potential are investigated in detail, in- 
cluding the distributions of densities and anomalous Hall 
currents. The quantized anomalous conductances appear 
in the band insulating regime at commensurate fillings. 
The magnitudes and chiralities of anomalous Hall cur- 
rents in the metallic regions cannot be described by the 
usual local-density-approximations. 

The rest of the paper is organized as follows. In Sec. 
HI! we give an introduction to the experimental setup and 
the orbital Zeeman coupling. In Sec. Mil a heuristic pic- 
ture is given to arrive at the Haldane model at large ro- 
tation angular velocities. In Sec. IIV1 the band structures 
including Berry curvatures and flat bands are studied. In 
Sec. El the spatial distributions of the particle density 
and anomalous Hall currents in the inhomogeneous har- 
monic trap is studied. Finite temperature effects are also 
studied. In Sec. I VII a brief discussion on the detection 
of the anomalous Hall current is presented. Conclusions 
are made in Sec. IVIII 



II. THE TIGHT-BINDING HAMILTONIAN 
WITH THE ON-SITE ROTATION 

In this section we describe the experimental setup by 
Gemelke et al. to realize the on-site rotation of opti- 
cal lattices [H-Q, and then construct the effective tight- 
binding model for such a system. 



A. The experiment setup by Gemelke et al. 

The honeycomb optical lattice was experimentally real- 
ized quite some time ago [13, [H, [zl, [73[ . It is constructed 
by three phase coherent coplanar laser beams with polar- 
ization along the z-axis, intersecting each other with 120° 
in the icy-plane. The schematic diagram of the experi- 
ment setup is shown in Fig. [TJ The potential minima of 
the interference pattern form the honeycomb lattice if the 
laser frequency is blue-detuned from the atom resonance 
frequency. The advantage of this technique is that the 
phase shift in the laser beams only leads to shift entire 
lattices without destroying the lattice geometry. 

The on-site rotation technique by Gemelke et al. was 
originally applied to the triangular lattice [H, Q . It would 
be straightforward to apply the same method to the hon- 
eycomb lattice. Two electro-optic modulators are placed 
in two of the laser beams, and the phase modulated po- 
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FIG. 1: (Color online) a) The sketch of the honeycomb optical 
lattice. The three laser beams cross each other at 120° in the 
icy-plane. Phase modulators are placed in the paths of the two 
beams, b) The illustration of the on-site rotation by Gemelke 
et al. 1-3]. The entire lattice takes the motions of a fast 
oscillation with the frequency wrf and a slow precession with 
the frequency Q, as schematically plotted with the red solid 
lines in one lattice site. After taking the time average of the 
fast oscillation, atoms feel that each site is rotating around its 
own center with the precession frequency Q, which is plotted 
with the red dot dash line around each site. 

tential is Q 

V(f) = Vb[cos(£i • f +<j> + ) + cos(fe 2 • r + (j>~) 

+ cos(fe 3 -r-0 + -0-)], (1) 

where fe, = h^ijk(Qj ~ 9fe)j <k{i — 1> 2, 3) are the wavevec- 
tors of the three coplanar laser beams satisfying \q\\ = 
1 92 1 = 1 93 1 =9- In the following of this paper, we use the 
definition of recoil energy E r = h 2 q 2 /(2M) where M is 
the atom mass. Please note that this definition of E r is 3 
times smaller than that used in Ref. [68[ which is defined 
as h 2 k 2 /(2M). In Eq. [TJ ^ = rjsm(nt± ^)sin(w RF £) 
where fl is the slow precession frequency; r\ is a phase 
modulating constant which determines the amplitude of 
the oscillation; wrf is the fast rotation frequency at radio 
frequency. Atoms do not follow the fast oscillation and 
only feel a time average of the potential as 

3 

V(r, t) = V J2 [Mt) cos(h ■ r)] , (2) 

i—l 

where A^[t) — J [i]sm(ilt + 3-)]; Jo is the zeroth order 
Bessel function; 77 is a small parameter Q. 



Eq. [2] still maintains the same lattice translational 
symmetry. Around each potential minimum r*o in the 
original lattice without rotation, the potential can be ex- 
panded to the second order, yielding a slightly anisotropic 
harmonic potential: 

V(r-r ,t) « _|_(l-_)| r _ nj | + — 

x |r-r | 2 cos[| + 2(m + ^_ Fo )]}, (3) 

which rotates with a slow frequency of fi. Here fr-r is 
the polar angle of r — f . The slight deformation of the 
optical potential processes around each site center, which 
can be regarded as an on-site rotation. 

B. The tight-binding Hamiltonian 

Now we construct the effective tight-binding model to 
describe the above system. First of all, each lattice site 
is rotating around its own center, and there is no overall 
rotation of the entire lattice. In other words, the sys- 
tem still has the lattice translational symmetry. There 
should be no vector potential for inter-site hopping asso- 
ciated with the Coriolis force. Within each site, the rota- 
tion angular velocity couples to the onsite orbital angular 
momentum through the orbital Zeeman coupling. Such a 
coupling also exists in solid state systems in the presence 
of external magnetic fields. However, the typical energy 
scales of the Zeeman couplings, including both spin and 
orbital channels, are at most at meV which are tiny com- 
pared to band widths. They usually do not change the 
band topology. The advantage of the experiments by 
Gemelke et al. [3-[3| is that the orbital Zeeman energy 
scale can easily reach the order of kHz, which is compa- 
rable to band widths. 

The orbital Zeeman term from the onsite rotation be- 
comes important in orbital bands with angular momenta 
higher than s. For the p^^-orbital bands, one of us [56[ 
introduced the following coupling as 

H L = -fi L z{r) = {pl,fPy,f ~ Pl,fPx,f} ■ ( 4 ) 

ft ^ 

It breaks the degeneracy between p x ± ip y states, and 
induces topologically non-trivial band structures as pre- 
sented in later sections. 

The remainin g pa rt of the tight-binding Hamiltonian is 
as usual. In Ref [66l.l68j. one of us studied the p^^-orbital 
bands in the honeycomb optical lattice filled with spin- 
less fermions, which is the counterpart of graphene de- 
scribed by the p z -orbital but exhibits fundamentally dif- 
ferent properties. The tight-binding Hamiltonian reads 

3 

H o = *|| Y {plfP*-f+e z + h.c.} - fi Y n (r)l 5 ) 

r£A i=l r£A(BB 

where ei : 2 = ±yi + \y and £3 = — y; A and B are 
indices of two different sublattices; t\\ is the c-bonding 
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FIG. 2: (Color online) The pattern of the induced complex- 
valued NNN hopping at Q ^> t« , which is generated by the vir- 
tual hop ping between orbitals with opposite chiralities. From 
Wu [H]. 

describing the longitudinal banding of p-orbitals along 
the bond direction; /i is the chemical potential; n{r) = 
P x pPx.f + Py pPy,r is the filling number at site r. The 
operators pi t ? are defined as the projection of p-orbital 
along the vector <§j as pp = p X; fX + p y .py. Rigorously 
speaking, t\i should be time-dependent which depends on 
the oscillation amplitude rj. Here we neglect this time- 
dependence by assuming rj is small, tn is positive as a 
result of the odd parity of the p -orbitals. The 7r-bonding 
t± is much weaker than the er-bonding. For example, 
t±/t\\ can be easily suppressed around 1% [56( within 
realistic experimental parameters of Vo/E r — 15, thus 
the t± is not considered in most of this paper unless in 
Sect. IIV CI The band Hamiltonian to be investigated 
below is the combination between Eq. |4] and Eq. [5] as 

H = H +H L . (6) 



which leads to the QAHE with v = ±1. In the system 
with the open boundary conduction, unidirectional edge 
currents appear surrounding the system, i.e., the edge 
currents are chiral. 

Before the detailed study on the band structures of 
our p-orbital Hamiltonian Eq. [51 we present an intuitive 
picture that the on-site rotation induces complex- valued 
NNN hopping terms in the limit of fi/iii ^> 1 as in the 
Haldane model. As a result, the Gemelke type rotation 
provides a possibility to realize the QAHE state in the 
cold atom experiments. In the presence of rotation, the 
on-site eigen-orbitals become p x ± ip y with an energy 
splitting of 20. When ft tn, each level of p x ± ip y 
broadens into a band without overlapping each other. 
We consider the case that fl > 0, such that the low en- 
ergy sector of the Hilbert space consists of the p x + ip y 
orbital state. The leading order term of the effective 
Hamiltonian in this sector is just the NN hopping with 
the hopping integral of ^iy. Moreover, the second order 
perturbation process generates the NNN hopping with 
complex-valued integral as explained below. 

Let us consider the two-step virtual hopping process 
illustrated in Fig. [5] In the first step, the atom starting 
from the low energy sector of the p x +ip y orbital in the A- 
site hops into the high energy sector of the p x —ip y orbital 
in the nearest neighbor B-site. The phases along the AB- 
bond is 30° from the A-site and 150° from the B-site, thus 
there is a phase mismatch of 120°. The corresponding 
hopping integral is complex-valued with \t\\e l 3 n . Sim- 
ilarly, during the second step the atom hops back into 
the p x + ip y orbital in the NNN A-site with the com- 
plex hopping integral |t||e l i 7r . The hopping process is 
(A+) — > (B-) — > (A'+) where ± represents the chiral- 
ity of p x ± ip y orbitals. The corresponding amplitude is 
calculated as follows: 



III. THE APPEARANCE OF THE HALDANE 
MODEL AT LARGE ROTATION ANGULAR 
VELOCITIES 

Haldane proposed a tight-bonding model for the 
QAHE effect whose Bloch wave band structure is topo- 
logically non-trivial [13] • The Hamiltonian of the Hal- 
dane model is defined in the honeycomb lattice, which 
reads 

H = -tJ2 {a\aj + h.c.} + {''.."-", + h-c], (7) 

where (ij) represents the nearest-neighbor (NN) hopping; 
((ij)) represents the next-nearest-neighbor (NNN) hop- 
ping. The NNN hopping t' — |£'|e ±l * is complex-valued 
whose argument takes ±</> if the hopping from i to j is 
anticlockwise (clockwise) with respect to the plaquette 
center. Eq. [7] breaks time- reversal symmetry. The band 
spectra exhibit two gapped Dirac cones in the Brillouin 
zone (BZ), whose mass values have opposite signs. The 
band structure has the non- vanishing Chern numbers ±1, 



(A' + \H \B-)(B-\Ho\A+) 

= 

= -1^ /3 - (8) 

All the NNN hoppings have the same phase value fol- 
lowing the arrows, which is exactly the same as in the 
Haldane model. The above analysis applies to the high 
energy sector as well. Thus we have two copies of the 
Haldane model, each for the p x ± ip y bands, respectively. 



IV. BAND STRUCTURES IN THE 
HOMOGENEOUS SYSTEM 

In this section, we present the band spectra in the 
homogeneous system with the periodical boundary con- 
ditions (PBC). The general structure is studied in Sec. 
IIV Al and the interesting flat band structure is presented 
in Sec. HVBl 
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(a) O /f,|=0 



(b) n/fn=l/2 



where H a b{k) is written as 




FIG. 3: (Color online) The band structures of the Hamilto- 
nian Eq. [IT]. With increasing rotation velocity from (a) 
n/tii = to (f) fl/tu — 1.7. The flat bands appear at 



Q/tu 



Energies are in unit of hopping strength tii and 



wavevectors are in the unit of 1/a, which is the reciprocal of 
the lattice constant. The energies of these two flat bands are 
E/t«=±l 



A. The general band structures 

We define the four-component spinor representation 
for the two-orbital wavefunctions in the two sublattices 
as 



\lp%) = PAx(k),PAy(k),PBx(k),PBy(k) 



(9) 



After performing the Fourier transform, the Hamiltonian 
Eq. [6] becomes 



(10) 



/ tQ 
-iQ. 
h.c. 
\ 



>/3 ^gifc-e*! 



o 

-iQ, 



iQ 
0, 



(.11) 



The band structure in the absence of rotation i.e., 
ft = 0, has been studied in Ref. [HI, HH, which includes 
both flat bands (the bottom and top bands) and two dis- 
persive bands with Dirac cones as depicted in Fig. [3] 
(a). The flat bands and dispersive bands touch at the 
center of the first BZ, and two dispersive bands touch 
at Dirac cones. The location of the Dirac cones are at 
k = (±3^,0) [gi, [68|. The band flatness means that 
the corresponding band eigenstates can be constructed as 
localized states in real space. Each hexagonal plaquette 
supports one localized eigenstate whose orbital configu- 
ration on each site is alon g th e tangent direction as pre- 
sented in Fig. 2a in Ref. [66|. When the filling is inside 
the flat bands, interaction effects are non-perturbative. 
This results in the exact solutions of the Wigner crys- 
tallization for spinless fermions [66]] and the flat band 
ferromagnetism for spinful fermions 69]. For the dis- 
persive bands, although their spectra are the same as 
in graphene, their eigen- wavefunctions are fundamentally 
different exhibiting rich orbital structures as presented in 
Ref. [H. 

When the on-site rotation is turned on, i.e., ft > 0, 
band gaps open. The previous touching points between 
the first and second bands at f2 = split. The lowest 
band is no longer flat, and the center of the second band 
is pushed up as depicted in Fig. [3] (b). The Dirac cones 
between the middle two dispersive bands also become 
gapped. In this case, the topology of each band n (n — 
1 ~ 4) is characterized by the Chern number defined as 



C»i — I d k F n ,xy{k) , 



where the Berry curvature F ntXy is defined as 



(12) 



(13) 



and A n ^([i = x,y) is the Berry gauge potential in mo- 
mentum space defined as 



An^ = i(lp n (k)\dk JV'n(fc)) • 



(14) 



The Chern number patterns at Q > have been cal- 
culated in Ref. [56[, and the distribution of the Berry 
curvatures F n , xy in the BZ is depicted in Fig. 2 of Ref. 
[56| . Below a critical value of the rotation angular veloc- 
ity ficA|| — f i the Chern number pattern reads 



d = -C 4 = 1; C 2 = -Ca = 0. 



(15) 



At Sl c A || — §, a single Dirac cone connecting the second 
and third bands shows up at k = (0,0), which triggers 
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a topological phase transition. Beyond f2 c , this Dirac 
point becomes gapped, and the Chern number pattern 
becomes 

Ci = -C 2 = 1; C 3 = -C 4 = -1. (16) 

In this case, the band structure is qualitatively the same 
as the two copies of the Haldane model as discussed in 

sec. nn 

B. Flat bands at Q./t\\ = § 

It is evident in Fig. [3] (c) that the second and third 
bands become flat at O/in = |. In this part, we discuss 
various properties of the flat bands, including the local- 
ized eigen-states, the distribution of the Berry curvature, 
and the interaction effects. 



1. Localized eigen-states 

The band flatness usually means that the eigen-states 
can be reconstructed as localized states in real space. 
We assume that each localized eigenstate exists within a 
single hexagon plaquette constructed as follows 

6 

fe) =E(-)^ le4 °" 1)0 { cos ^l^) "Sin 

(17) 

where R is the coordinate of the plaquette center; j is the 
site index within the same plaquette and d 3 ■ = (j — l)-f ; 
e 1 ^ is the phase factor to be determined satisfying the 
periodical boundary condition e l6< ^ = 1; the factor of 



(— 1 is a sign convention because of the odd parity of 
the p-orbitals. The p-orbital configuration on each site j 
is along the tangent direction. Substituting Eq. [IT] into 
the band Hamiltonian, we arrive at the condition for Eq. 
fTTl to be the eigenstate as 

Q = — —sine/), E — —-cos4>, (18) 

where 4> — 0, ±| , ±^-, 7r. For the cases of <f) = and 
7r, they are the situations studied before in Ref. (6(| 
without the on-site rotation. The other four cases are 
with the on-site rotation. Without loss of any generality, 
we take </> = — -| and 4> — — if sucn that ^A|| = § > 
and E/t\\ = ±|. The schematic diagram of these two 
typical localized state is shown in Fig. @] (a) and (b), 
respectively. 

The main difference between these two groups of local- 
ized states at 57 = and f2/i|| = | is that there exists 
a current around each plaquette for the latter case. The 
current operator along each bond is defined as 

Jr,r+e z = ^ ' $i) (Pf+ei ' - h-C.^. (19) 

For the localized plaquette eigen-states of both bands 
with E — the currents have the same value and 

chirality as 

J = -V3CI. (20) 

Eq. [20] indicates that the current direction is opposite 
to the rotation and the magnitude is proportional to the 
angular velocity f2. 

At = f j we solve the eigenvectors for the two fiat 
bands in momentum space as 



1 — ifc-ei=pixi 



N T (k) 



1 „-ifc-e 2 =BX2 
2 e 



V3 g-'ifc-ei=Fixi _ V_l c -'ifc-e?T'tY2 



V 



zp^v^3 g^-eiT^Xl _ V3gife-e 2 =F 



'X2 ~1 



(21) 



where ^>=p(fc) represent eigenvectors for the bands with 
E/t\\ = =f|, respectively; xi = |, X2 = ^, and %3 = if- 
The normalization factors N^z(k) read 



N^(k) = 2 



cos(fc • hi =F gTr) 



(22) 



These flat band Bloch wave states can be represented as 
the linear superpositions of the localized eigenstates in 
Eq. [T7]as: 

\^ k J = ^L=Y,^ k - S ^ R ) (23) 



N T (k) R 



where \ipR) is defined in Eq. [T71 



2. Brief discussions on interaction effects 

The band flatness means that interaction effects are al- 
ways important compared to the vanishing kinetic energy 
scale. In our previous studies [66], [68|, [6i| , we have exam- 
ined the non-perturbative effects in the flat bands (the 
lowest and highest bands) in the same system at fi = 0. 
In Ref. (6(1 [68[ , we have shown that the flat bands re- 
sult in the exact solution of Wigner crystal configuration 
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FIG. 4: (Color online) Configurations of the localized eigen- 
states of the two flat bands at Q/t\\ = § for (a) E/t\\ = — § 
and (b) E/t\\ = §, respectively, ui — e' s is the relative phase 
factor between neighboring sites. The plaquette currents di- 
rections are clockwise for both (a) and (b) in the opposite 
direction of the Q. 



(a) (b) 




FIG. 5: (Color online) The Berry curvature F^in unit of 
a 2 ) distribution at f2/tii ~ I f° r * ne ls t ( a ) an d 2nd (b) bands. 
The Chern numbers for the 1st and 2nd bands are 1 and 0, 
respectively. The unit of wave vector is 1 /a. 



3. Berry curvatures v.s. local eigen-states 

The current carried by the localized eigenstates of the 
flat bands depicted in Fig. @]is chiral. It looks very sim- 
ilar to the classic picture of cyclotron orbit of electrons 
in the external magnetic fields. We would expect that 
in a system with the open boundary condition, the fully 
filled flat band would result in edge currents and con- 
tributes to the quantized anomalous Hall conductance. 
However, we need to be very careful with this, which 
turns out to be incorrect. We have performed a pre- 
liminary diagonalization for a finite size system with the 
open boundary condition. The number of degeneracy for 
the flat bands equals to the number of plaquettes plus 
1. We conjecture that this extra state should not belong 
to a particular plaquette but rather distribute along the 
edge, which carries a current in the opposite direction 
and cancels the contribution from other plaquette states. 
Further examinations on this problem will be deferred to 
a later publication. 

We calculate the Berry curvature distributions at 
f2/i|| = I f° r t nc 1st and 2nd bands as presented in 
Fig. [5] (a) and (b), respectively. Those of the 3rd (4th) 
band are just with an opposite sign compared to the 2nd 
(1st) band due to the particle-hole symmetry of the band 
Hamiltonian. The 1st and 4th bands are topologically 
non-trivial with the Chern number ±1. However, the 
2nd and 3rd bands, which are flat, are topologically triv- 
ial with the zero Chern number. In fact, these two bands 
should not contribute to quantum anomalous Hall con- 
ductance when they are fully filled. This is confirmed 
from the anomalous Hall current calculation in the inho- 
mogeneous trap as presented in Sec. |V] 



for spinless fermions in the lowest band. At (n) = g, 
which corresponds to that | of the fiat band plaquette 
states are occupied, the occupied plaquettes form a trian- 
gular lattice structure without touching each other. As 
filling increases, exact solutions are no longer available. 
Self-consistent mean-field theory calculation shows a se- 
rials of insulating states with different orbital orderings 
at commensurate fillings 



Similarly, in Ref. [6° 



we found the exact flat-band ferromagnetism for spinful 
fermions in the flat bands. 

For the flat bands occurring at il/t^ = |, the physics 
will be similar to the previous studies at tl — 0. How- 
ever, the flat bands here are in the middle. When /i lies 
in the flat bands, there are always background particles 
or holes filling in the dispersive bands. The solutions 
of the Wigner crystal and flat-band ferromagnetism is 
only valid if the interaction energy scale is smaller than 
the band gaps between the flat bands and the dispersive 
bands. For example, when the filling is inside the second 
band, the effect from the background filling cannot be 
neglected, if the interaction energy scale is stronger than 
the band gap. 



C. Effects of the 7r-bonding to band structures 

So far, we have neglected 7r-bonding t± which can be 
easily suppressed around 1% of tu at intermediate optical 
potential strength (68|. Here we explicitly present its ef- 
fects to band structures for the case of a realtively weak 
lattice potentials by taking t±/t\\ = 0.05. The projec- 
tions of p^y-orbitals perpendicular to the ei.2,3 directions 

are defined as p' 12 — ~\Vx ± g Pti'Ph = respectively. 
The 7r-bonding Hamiltonian can be written as 



E {p'liPWa^i + h-c), (24) 



where the hopping integral of the 7r-bonding has the 
opposite sign to that of the cr-bonding. In momentum 
space, Eq. [M] transforms into 

H v = -t ± Y,^ a (k)H^ a p{k)^p(k), (25) 
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FIG. 6: (Color online) The band structures with both tt- 
bonding (t±/tn — 0.05) and on-site rotation Q. Only the 
lower two bands are presented, and the spectra of the other 
two bands are symmetric with respect to zero energy, (a) 
At fl/t|| = 0, the two bands remain touching at the center 
of the BZ. The lowest band develops dispersion at the order 
of iii. (b) At fl/tu = 0.2, the gap opens between the lower 
two bands. The lowest band is topologically non-trivial and 
nearly flat. 



with the martix kernel H n (k) as 



( §(e iS ' 41 +e i&82 ) + , 

h.c. 



,ik-e 3 3^_ e ifc-e! 








+ e 1 ' 



) 



(26) 



The effects of 7r-bonding t± are presented in Fig. [6] 
The spectra remain symmetric with respect to zero en- 
ergy, and thus only the lower two bands are presented. 
As presented in Ref. (68[, at f2 = 0, the bottom bands 
are no longer rigorously flat but develops a finite width 
at the order of t±. The lower two bands remain touch- 
ing at k = (0, 0) with parabolic spectra, and the bottom 
band has a negative curvature. With increasing f2, as in 
the case of t± = 0, the band gap at the order of 51 opens. 
Furthermore, f2 lowers the energies of the bottom band 
near the center of the BZ, which suppresses its dispersion. 
As a result, we arrive at a nearly fiat band with non-zero 
Chcrn number. The ratio between the width of the bot- 
tom band and the gap between the lower two bands can 
reach the order of 5 as shown in Fig. [5] b. Recently, we 
notice that the nearly flat bands with non-trivial Chern 
number have been attracting attention, for its possible 
realization of fraction quantum Hall states in the lattice 

MM. 



ANOMALOUS HALL CURRENTS IN 
HARMONIC TRAP POTENTIALS 



In section IIV1 the homogeneous p-orbital system with 
the PBC has been studied, in which the wavevector k is 
a good quantum number. However, in reality the honey- 
comb lattice is inhomogeneous with a soft harmonic con- 



LU 




0.03 0.06 0.09 

g(E) 

FIG. 7: (Color online) The DOS of the p^-orbital model 
Eq. [5] and Eq. [4] hi the homogeneous system. The red solid 
lines are the DOS at £l/t\\ = § whereas the blue dash lines 
are the DOS at n/ty = §. 



fining trap. In this section, we shall consider the anoma- 
lous Hall currents in such a realistic system. 

The trapping potential adds a new term in the band 
Hamiltonian Hq + Hl of Eqs. [5] and H] as 



H T = ^V T {r)n{r). 



The trapping potential Vr(r) reads 



(27) 



V T (r) = -Muj^r 2 
where a is the lattice constant, 



2 V ' 



(28) 



in 



" 



Mujt is the trapping length scale. The typical value 

of the trapping frequency lot is in the order of 10Hz, 
and that of the recoil energy E r is roughly several kHz 
77|. In Ref. [Hj], we have calculated that t\\/E r = 0.24 
for V~o/E r = 15, thus hu/t\\ is at the order of 0.1. The 
typical trapping length scale is several lattice constants. 
Taking into account all these factors, we choose a conve- 
nient value of [5 = 0.01 for later calculations. 

In the inhomogeneous system with trapping potential, 
the on-site rotation induces the circulating currents along 
the azimuthal direction. We will study the spatial dis- 
tributions of the these anomalous Hall currents and par- 
ticle density. Because the band topology has a transi- 
tion at fi c /t|| = |, the results are presented at different 
sets of parameters below, at, and above Q c . Our results 
are calculated by using the eigen-wavefunctions from the 
numerical diagonalization of the free Hamiltonian in an 
open lattice with the trapping potential. We also use 
a modified local-density-approximation (LDA) to under- 
stand the exact results. The size of the lattice is with the 
radius of r/a = 40. 
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■ : (Color online) The particle density distributions 
v.s. the radius r at (a) Q/ty = | and (b) §, which 
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exhibit a four-layer wedding cake shape (the blue dots). The 
radius is in the unit of the lattice constant a. The DOS at 
the local chemical potential ni oc {r) in the LDA approxima- 
tion is plotted with the black dash lines. /j,i oc (r)/tu — 2.3 at 
the center of the trap and uir/tn = 0.1. 
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FIG. 9: (Color online) The effective local Hall conductance 
&po(r) v.s. r defined in Eq. [31] at Q/t\\ = ^ (a) and § (b). 
The radius is in the unit of the lattice constant a. The results 
from diagonalizing the free Hamiltonian with the trapping po- 
tential is marked with asterisks, and those from the modified 
LDA are plotted with dashed lines. o p g is quantized in the 
insulating plateaus with commensurate fillings of (n(r)). 



A. Low rotation angular velocity 

In this subsection the angular velocities are taken as 
f2/t|| = \ and | below 51 c . The chemical potential is cho- 
sen as = 2.3 to guarantee that all bands are filled 
at the center of the trap. The spatial distribution of the 
particle density exhibits a four-layered wedding cake-like 
structure. The density plateaus correspond to the band 
insulating regions, where the local chemical potential, 
defined as fJ,i oc (r) = /i — Vr(r), lies inside band gaps. 
Furthermore, the anomalous Hall currents flow along the 
tangent direction, whose conductances are quantized in 
the insulating regions. 



1. The insulating plateaus of (n(r)) 

We first present the density of states (DOS) for the 
p-orbital bands in the homogeneous system at fl/t^ = -|, 
and § in Fig. [7] The DOS is defined as 

9(E) = J -^5{E{k)~E). (29) 

At ri/t|| = j, the strong divergence of the DOS is indi- 
cated at the 2nd and 3rd bands due to the appearance of 
the flat bands. At the small value of O/tii = i, the DOS 
of the 1st and 4th bands are larger than the 2nd and 3rd 
bands, which is a reminiscence of the band flatness at 
n = 0. Moreover, it is obvious that the band gaps open 



at tt c > 0. For the chemical potential fi lying in the band 
gaps, the system is in the band insulating states with 
the commensurate values of the particle number per site 
(ri) = |, 1, |, and 2, respectively. 

In the inhomogeneous trap, the real space distribu- 
tions of the filling number (n(r)} are calculated by us- 
ing the eigenstate wavefunction obtained through diago- 
nalizing the Hamiltonian, which are depicted in Fig. [8] 
for = i and |. In both cases, plateaus appear at 

(n) = |, 1, |. These plateaus can be understood within 
the LDA picture. Recall the band structure in Fig. [3]and 
the DOS in Fig. [7] When /x; oc (r) lies in the band gaps, 
the filling number stops increasing until /i/ c(?') reaches 
next band edge. The local DOS at site r at the en- 
ergy jJ,io C {r) is also plotted, which roughly proportional 
to d r n(r). It is clear that the locations of the plateaus of 
(n(r)) are and band gaps are consistent. 

Because the honeycomb lattice breaks the SO(2) ro- 
tational symmetry down to the 6-fold one, lattice sites 
with the same magnitude of r may have different values 
of (n(r)). They are slightly scattered in the metallic re- 
gions between different plateaus as depicted in Fig. [8] 
(a) and (b). In the case of n/ty = |, the distribution 
of (n(r)) exhibit devil's stair-like features [78|, [79J. In 
the cliffs as filling in the middle two flat bands whose 
degeneracies are slightly lifted by the potential gradient. 
If with interactions, the flat band regime may further 
exhibit plateaus of Mott-insulating states with orbital 
orderings, which will be deferred to a later research. 
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FIG. 10: (Color online) The pattern of the anomalous Hall 
current in the honeycomb lattice with the confining trap at the 
rotation of fl/tii = 0.75. The blue dotted (red solid) lines rep- 
resent the counterclockwise (clockwise) anomalous Hall cur- 
rents, respectively. The color depth indicates the magnitude 
of the current. The reversed direction of the anomalous Hall 
currents (red solid lines) in the metallic regions between two 
neighboring plateaus can be explained as the anomalous con- 
tribution from the gradient of (n(r)). 

2. QAHE currents in the insulating density plateaus 

Due to the non-trivial topology of the band structure, 
anomalous Hall currents circulate along the azimuthal di- 
rection due to the radial potential gradient. The plateaus 
of the filling number (n(r)) correspond to the insulat- 
ing quantum anomalous Hall regions with quantized Hall 
conductance. Compared with the usual quantum Hall 
systems, the on-site rotation breaks time-reversal sym- 
metry and brings non-trivial topology to the band struc- 
tures without Landau levels. 

The anomalous Hall current along each bond is calcu- 
lated by using the eigen-wavefunctions obtained from the 
diagonalization of the real space p-orbital Hamiltonian as 

J?,r+&i = i-r^iMPp- g i)(Pf+^ ■ e») - h.c.\), (30) 

n 

where (| |) represents the ground state at T = or thermal 
average at finite temperatures. Let us focus on those 
bonds orienting along the azimuthal direction, and define 
the effective local Hall conductance as 



where r is the radius of the middle point of the bond; 
jo is the current density; Vt is the trapping potential. 
In our honeycomb lattice system, the current density is 
defined as the current on each bond Jfft+^ ( divides the 
distance between neighboring parallel bonds. 



In the homogeneous systems, the Hall conductance is 
represented as 

a pe = ]l2^22 d 2 kF ltXV (k)n f (i,k), (32) 

i 

where n/ is the Fermi distribution function; i is the band 
index. When the chemical potential is inside band gaps, 
o pQ is quantized as the sum of the Chern numbers of the 
occupied bands [H [H Q 

°p» = \X' Ci ' (33) 

i 

For the cases of = \, |, the Chern number patt ern is 
the same as C\ = — C 4 = 1 and C 2 = — C 3 = 56]. The 
quantized Hall conductances reads 0, 1, 1 and 1 as \x lies 
from above the band top down to the three consecutive 
three band gaps. 

The results of cr°g(r) (defined in Eq. |3"TT) v.s r are 
marked as asterisks in Fig. [S] (a) and (b) for f2/i|| = \ 
and §, respectively, which are obtained by diagonaliz- 
ing the free but inhomogeneous Hamiltonian. The real 
space circulating current pattern at fi/iy = f is depicted 
in Fig. [TO] The quantized Hall conductances in the in- 
sulating plateau regions can be understood within the 
LDA picture. At the center, the local chemical poten- 
tial ni oc (r) lies above the band top, and the conductance 
is therefore zero. As moving into the insulating density 
plateaus of (n(r)) = |,1, and ^, a^g is quantized at 
1/h. Counterclockwise currents are plotted as blue dash 
line under the harmonic trap potential. When the radius 
r > 30, the Fermi level is lower than the band bottom, 
thus the current vanishes again. 



3. Anomalous Hall currents in the metallic regions 

Between two adjacent insulating plateaus, the system 
is metallic with incommensurate fillings of (n(r)), there- 
fore the anomalous Hall conductances are non-quantized. 
The Hall current response to the radial potential gradi- 
ent is non-local in these inhomogeneous metallic regions. 
The effective Hall conductance a p g defined in Eq. [31] can- 
not be obtained from Eq[35]of the homogeneous system 
by using LDA with a local chemical potential /z; oc (r). 
For example, the currents can reverse the direction to 
be clockwise in the metallic regions in Fig. [H] (a) and 
(b). However, because of the Chern number pattern of 
C2 = C3 = and C\ = — C4 = 1, the Hall conductance 
a p g defined in Eq. [32] is always positive within the LDA 
at ~ < (n) < 1 which corresponds to that /i; oc (r) lies in 
the gap between the 1st and 2nd bands, The naive LDA 
results would only give rise to counterclockwise Hall cur- 
rents in these two metallic regions. 

Now we propose a modified LDA method to fit the 
above exact results from diagonalization. We define the 
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FIG. 11: (Color online) The particle density distributions 
(n(r)} v.s. the radius r at (a) Q./t\\ — | and (b) 2 (the 
blue dots). The radius is in the unit of the lattice con- 
stant a. The DOS at the local chemical potential fj,i oc (r) in 
the LDA approximation is plotted with the black dash lines. 
MiocM A|| = 3-5 at the center of the trap and uJT/t\\ = 0.1. 
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FIG. 12: (Color online) The effective local Hall conductance 
o-po(r) v.s. r defined in Eq. [31] at £l/t\\ = § (a) and 2 (b). 
The radius is in the unit of the lattice constant a. The results 
from diagonalizing the free Hamiltonian with the trapping po- 
tential is marked with asterisks, and those from the modified 
LDA are plotted with dashed lines. 



effective driving force as the derivative of the spatial- 
dependent part of the ground state energy density as 



F(r) 



1 d 



{(> ioc (r) - £ B (Q))n(r)} 



n(r) d r 

Adrift + F diS , (34) 



where Eb(Q) is the band bottom energy; Firift and -Fdiff 
are defined as 



-Fdrift = -d r V T (r), 



Fan 



n(r) 



d r n{r). 



(35) 



Fdrift comes from the gradient of the trapping potential, 
while Fdiff is the chemical pressure from the particle den- 
sity gradient. 

Correspondingly, the anomalous Hall currents can be 
interpreted by two contributions from the "drift" and 
"diffusive" Hall currents as 

Je(r) = Jdrift,fl(V) + JdiK,e(r) 

= cr p g(r){F dlift (r) + F di s(r)}, (36) 

where a p e (r) is obtained from Eq. [32] in the LDA by 
using the local chemical potential jLt; oc (r). Thus u c p ^(r) 
defined in Eq. [31] is related to & P e(r) through 



eff / \ , s f, , Fdiff(r) -I 



(37) 



The results of the effective Hall conductance cr°* (r) using 
this modified LDA is presented in Fig. [H] with dashed 
lines, which nicely agrees with the exact results. 



In the insulating plateaus, .Fdiff = 0, thus the modi- 
fied LDA reduces back to the naive LDA. However, in 
the metallic regions, d r n(r) has the opposite direction to 
d r Vr- As a result, the direction of JdiS,e is also oppo- 
site to that of Jdrift,e- The reversed direction of the Hall 
currents in the metallic regions can be understood as the 
contribution of "diffusive" Hall current dominates over 
that of the "drift" Hall current. 



B. Large rotation angular velocities (fl/iii 



1,2) 



In this subsection, we consider the large rotation angu- 
lar velocities at f2/t|| — f an d 2 which are at and above 
fi c /tii = §, respectively. The band structure topology 
above f2 c changes to a different Chern number pattern 
of Ci = — C2 = C3 = — C4 = 1. The chemical potential 
n/t\\ is chosen to 3.5, which guarantees that all bands are 
filled at r = 0. 

The distributions of the filling number (n(r)) are de- 
picted in Fig. E] at (a) fi/t|| = § and (b) ft/ty = 2, 
respectively. At Sl/i|| = §, the 2nd and 3rd bands touch 
each other at a Dirac cone located at the center of the BZ. 
The DOS vanishes linearly, and thus the density profile 
exhibits a soft slope instead of a flat plateau in Fig. [TT] 
(a). In both Fig. [TT](a) and (b), the density distributions 
between the 1st and 2nd bands also exhibit soft slopes 
although there do exist a band gap in the homogeneous 
system. This is because the potential gradient increases 
as r goes larger in the confining trap, which closes the 
small gap between the 1st and 2nd bands at large values 
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FIG. 13: (Color online) The radial distributions of (a) the 
anomalous Hall conductance a p e(r,T) and (b) the entropy 
S(f,T) with Q/t\\ = f. The black dot, blue "+", magenta 
diamonds and red asterisks of the data represent temperatures 
at T = 0, T/t\\ = 0.1, T/t\\ = 0.2 and T/t\\ = 0.5, respectively. 



of O. 

The local anomalous Hall conductances defined in Eq. 
[31] are depicted in Fig. Q2]at f2/£|| = § an d 2. In the in- 
sulating plateaus between the 3rd and 4th bands, cr p g is 
close to the quantized value of 1 /h for both rotation angu- 
lar velocities. The small deviation comes from the finite 
width of the insulating regions. In the region with soft 
slopes of the distributions of (n(r)) between the 1st and 
2nd bands, a p e is significantly smaller than 1/h because 
this region is not rigorously insulating. In the inhomoge- 
neous metallic regions, the values of the local anomalous 
Hall currents are non-quantized which are determined by 
the combined effects from the gradients of the trapping 
potential and the density distribution. 



C. Temperature effects 

In this subsection, we briefly discuss the finite temper- 
ature effects to the anomalous Hall conductance. The 
QAHE is a topological property existing in band insulat- 
ing regions, thus it is robust against finite temperatures 
provided their scale is small compared to the band gap. 
On the other hand, we do expect that the anomalous Hall 
conductance in the metallic regions will be significantly 
affected by finite temperatures. 

The radial distribution of the local anomalous Hall 
conductance cr°f (r) v. s. r is plotted in Fig. [T3] (a) at 
different temperatures. The rotation angular velocity is 
taken as il/tn — §. In this case, the band gaps are at the 
same order of tn as shown in Fig. [7] a p e remains nearly 
quantized in the insulating regions for T/tu = 0.1 ~ 0.2. 



According to the calculation of band structures in Ref. 
68], i|| « 0.24E r at V /E r = 15 and the typical energy 
scale of E r is 0.1 ~ 0.2/iA", thus the QAHE signature 
should survive at the order of lOnK, which is an exper- 
imentally accessible temperature scale. In the metallic 
regions, naturally <x p g's are more strongly affected by fi- 
nite temperatures. As T further increases to the half 
value of tii , the quantized signatures of a p g disappear. 

We also present the entropy distributions in real space 
at various temperatures as shown in Fig. [T3] (b). The 
local entropy is defined as 

s(f,T) = -k B i iMoi a {»*cn ln n i( T ) 

i 

+ (l-riiCrjjHi-niCr)]}, (38) 

where kg is the Boltzmann constant; the subscript i 
is the index of energy levels; ipi(r) is the wavefunction 
at the location r; n^(T) is the Fermi distribution func- 
tion. At low temperatures (e.g. T/t\\ =0.1 and 0.2), 
S(f, T) concentrates in the gapless metallic regions, and 
remains negligible in the band insulating regions. We 
have the coexistence of the insulating QAHE regions and 
the metallic regions to hold a significant amount of en- 
tropy. The tolerance of the large residue entropy den- 
sities in the trap greatly facilitates the experimental re- 
alization of the QAHE state. As temperatures go high, 
say, T/tu = 0.5, the entropy distribution becomes more 
uniform, and there are no clear distinctions between in- 
sulating and metallic regions any more. This agrees with 
the picture in Fig. [T3] (a), in which the plateaus of quan- 
tized anomalous Hall conductance disappear. 



VI. EXPERIMENTAL DETECTIONS 

In experiments, the plateaus of commensurate fillings 
of atoms in Fig. [5] and Fig. [TT]can be observed by mea- 
suring the in-trap density distribution or the compress- 
ibility of the lattice system, as clearly demonstrated in re- 
cent experiments |80l - [83l |. However, the density plateaus 
cannot distinguish the conventional band insulators and 
the quantum anomalous Hall insulators. 

In solid state systems, the Hall conductivity are ob- 
tained from transport measurements, which are very dif- 
ficult for the cold atom experiments. Nevertheless, it 
has been proposed to detect through the response of the 
atom density to an external magnetic field (84|, which 
can be realized by further rotating the harmonic trap 
[8^ or coupling atoms with additional laser fields. In 
particular, the motion of atoms in laser fields leads to 
an artificial magnetic field, which has been observed in 
a recent experiment [86j. In the presence of an artificial 
magnetic field, the quantized anomalous Hall conductiv- 
ity cr xy = (^5) , according to the well-known Streda 
formula [13] derived for the quantum Hall effects in the 
solid state. Therefore the density of atoms changes lin- 
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early with respect to the applied magnetic field when a xy 
is quantized in some regions of the harmonic trap. 

Another possible method to detect the anomalous Hall 
current is as follows. We assume all atoms are initially 
prepared in a hyperfine ground state |1). To detect the 
anomalous Hall current shown in Fig. 9, we apply a local 
two-photon Raman transition using two co-propagating 
focused laser beams in a small area S to transfer atoms in 
S to another hyperfine state |2). A subsequent time-of- 
flight measurement of the velocity distribution of atoms 
in state |2) gives the initial velocity distribution (thus the 
current) of atoms in the state |1) in the optical lattice. 
The above density and current measurements provide the 
experimental signature of the quantum anomalous Hall 
effects in the p-orbital honeycomb lattice. 

VII. CONCLUSIONS AND OUTLOOK 

In summary, we have proposed the realization of the 
quantum anomalous Hall states in the cold atom optical 
lattices based on the experimentally available technique 
of the on-site rotation developed by Gemelke et al. This 
rotation generates the orbital Zeeman coupling whose en- 
ergy scale can reach the order of the band width. In the 
p-orbital bands of the honeycomb lattice, the band struc- 
tures become topologically non-trivial at any nonzero ro- 
tation angular velocities. A topological transition occurs 
at fi c /tii = | with different band Chern number pat- 
terns below and above f2 c . At O > Cl c , the band topol- 
ogy is equivalent to a double copy of Haldane's quantum 
anomalous Hall model. Flat band structures are also 
found at fi/iy = § whose localized eigenstates can be 
constructed as circulating plaquette current states. The 
flat band structures may bring strong correlation effects, 
such as Wigner crystal and ferromagnetism, when inter- 
actions are turned on. 



The effects of the spatial inhomogeneity to the p- 
orbital quantum anomalous Hall states are also investi- 
gated. At each commensurate filling of |,1,|, and 2, 
the density profile exhibits insulating plateaus, whose 
Hall conductances are quantized at integer values. In 
the metallic regions between two adjacent plateaus, the 
anomalous Hall currents are determined by the non-local 
response, which can be understood as the combined ef- 
fects of the gradients of the confining potential and par- 
ticle density. We have also showed that the QAHE is 
robust at finite but low temperatures compared to band 
gaps. 

We further point out that the generation of the quan- 
tum anomalous Hall states from this "orbital Zeeman" 
effect is very general, not just for the honeycomb lat- 
tice. The advantage of the p-orbital honeycomb lattice is 
that an infinitesimal value of fi is enough to generate the 
quantum anomalous Hall states. For other generic lattice 
structures, beyond a critical value of f2 which is compa- 
rable to the band width, the orbital Zeeman effect gen- 
erates inverted orbital bands of different orbital angular 
momenta. The further hybridization among them brings 
non-trivial band topology, which is a similar mechanism 
to achieve topological insulators in semi-conducting sys- 
tems through spin-orbit couplings. A systematic study 
will be presented in a later publication. 

Acknowledgement 

M. Z. acknowledges Prof. Shi-qun Li for the sup- 
port and thanks Wei-Cheng Lee for helpful discussions. 
H. H. H. and C. W. are supported by NSF under No. 
DMR- 0804775 and AFOSR-YIP, M. Z. is supported by 
the NFRP-China Grant(973Project) Nos.2006CB921404, 
and C.Z. is supported by ARO (W911NF-09-1-0248) and 
DARPA-YFA (N66001-10-1-4025). 



[1] N. Gemelke, E. Sarajlic, and S. Chu, Rotating Few-body 
Atomic Systems in the Fractional Quantum Hall Regime, 
arXiv: 1007.2677 (2010). 

[2] N. Gemelke, Ph.D. thesis (2007). 

[3] E. Sarajic, N. Gemelke, S.-W. Chiow, S. Herrman, 
H. Miiller, and S. Chu, Coherent control of ultracold mat- 
ter: Fractional quantum hall physics and large-area atom 
interferometry (2009). 

[4] R. Karplus and J. M. Luttinger, Phys. Rev. 95, 1154 
(1954). 

[5] J. Smit, Physica 24, 39 (1954). 

[6] L. Berger, Phys. Rev. B 2, 4559 (1970). 

[7] T. Jungwirth, Q. Niu, and A. MacDonald, Phys. Rev. 

Lett. 88, 207208 (2002). 
[8] N. Nagaosa, J. Phys. Soc. Jpn 75, 42001 (2006). 
[9] N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 

N. P. Ong, Reviews of Modern Physics 82, 1539 (2010). 
[10] D. Xiao, M.-C. Chang, and Q. Niu, arXiv:0907:2021 

(2009). 



[11] Z. Qiao, S. A. Yang, W. Feng, W. Tse, J. Ding, 
Y. Yao, J. Wang, and Q. Niu, Chern Number Cre- 
ation in Graphene from Rashba and Exchange Effects, 
arXiv:1005.1672 (2010). 

[12] T. Tomizawa and H. Kontani, Phys. Rev. B 80, 100401 
(2009). 

[13] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 
M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

[14] M. Kohmoto, Ann. Phys. 160, 296 (1985). 

[15] R. Jackiw, Phys. Rev. D 29, 2375 (1984). 

[16] E. Fradkin, E. Dagotto, and D. Boyanovsky, Phys. Rev. 
Lett. 57, 2967 (1986). 

[17] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[18] M. I. Dyakonov and V. I. Perel, Physics Letters A 35, 
459 (1971). 

[19] J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

[20] J. Sinova, D. Culccr, Q. Niu, N. A. Sinitsyn, T. Jung- 
wirth, and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 
(2004). 



14 



S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, [52 
1348 (2003). 

Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. [53 
Awschalom, Science 306, 1910 (2004). 

J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, [54 
Phys. Rev. Lett. 94, 047204 (2005). 

B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science [55 
314, 1757 (2006). 

X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008). [56 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 [57 
(2005). 

D. N. Sheng, Z. Y. Weng, L. Sheng, and F. D. M. Hal- [58 
dane, Phys. Rev. Lett. 97, 036808 (2006). 

J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 
(2007). 

R. Roy, Phys. Rev. B 79, 195321 (2009). [59 
L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, [60 
106803 (2007). 

L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). [61 
H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. [62 
Zhang, Nat Phys 5, 438 (2009). 
C. Wu, B. A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. [63 
96, 106401 (2006). [64 

C. Xu and J. E. Moore, Phys. Rev. B 73, 064417 (2006). 
M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buh- [65 
mann, L. W. Molcnkamp, X. L. Qi, and S. C. Zhang, 
Science 318, 766 (2007). [66 

D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. 
Cava, and M. Z. Hasan, Nature 452, 970 (2008). [67 
D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, [68 
J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, 

et al., Science 323, 919 (2009). [69 
Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, 
A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, et al., Nat [70 
Phys 5, 398 (2009). 

Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. [71 
Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, 
et al., Science 325, 178 (2009). [72 
P. Roushan, J. Seo, C. V. Parker, Y. S. Hor, D. Hsieh, 
D. Qian, A. Richardella, M. Z. Hasan, R. J. Cava, and [73 
A. Yazdani, Nature 460, 1106 (2009). 

Z. Alpichshev, J. G. Analytis, J.-H. Chu, I. R. Fisher, [74 
Y. L. Chen, Z. X. Shen, A. Fang, and A. Kapitulnik, 
Phys. Rev. Lett. 104, 016401 (2010). [75 
T. Zhang, P. Cheng, X. Chen, J.-F. Jia, X. Ma, K. He, 
L. Wang, H. Zhang, X. Dai, Z. Fang, et al., Phys. Rev. [76 
Lett. 103, 266803 (2009). 
X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Phys. Rev. B 74, [77 
085308 (2006). [78 
C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, [79 
Phys. Rev. Lett. 101, 146802 (2008). [80 
M. Onoda and N. Nagaosa, Phys. Rev. Lett. 90, 206601 
(2003). [81 
R. Yu, W. Zhang, H.-J. Zhang, S.-C. Zhang, X. Dai, and 
Z. Fang, Science 329, 61 (2010). [82 
T. L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 
V. W. Scarola and S. D. Sarma, Phys. Rev. Lett. 98, 
210403 (2007). [83 
R. O. Umucalilar, H. Zhai, and M. O. Oktel, Phys. Rev. 
Lett. 100, 070402 (2008). [84 
S.-L. Zhu, H. Fu, C. J. Wu, S. C. Zhang, and L. M. Duan, 
Phys. Rev. Lett. 97, 240401 (2006). [85 
C. Zhang, arXiv:1004.4231 (2010). 



L. B. Shao, S.-L. Zhu, L. Sheng, D. Y. Xing, and Z. D. 
Wang, Phys. Rev. Lett. 101, 246810 (2008). 
T. D. Stanescu, V. Galitski, and S. Das Sarma, Phys. 
Rev. A 82, 013608 (2010). 

X.-J. Liu, X. Liu, C. Wu, and J. Sinova, Phys. Rev. A 
81, 033622 (2010). 

N. Goldman, I. Satija, P. Nikolic, A. Bermudez, M. A. 

Martin-Delgado, M. Lewenstein, and I. B. Spielman, 

Phys. Rev. Lett. 105, 255302 (2010). 

C. Wu, Phys. Rev. Lett. 101, 186807 (2008). 

G. Grynberg, B. Lounis, P. Verkerk, J. Y. Courtois, and 

C. Salomon, Phys. Rev. Lett. 70, 2249 (1993). 

P. Soltan-Panahi, J. Struck, P. Hauke, A. Bick, 

W. Plenkers, G. Meineke, C. Becker, P. Windpassinger, 

M. Lewenstein, and K. Sengstock, arXiv: 1005. 1276 

(2010). 

W. V. Liu and C. Wu, Phys. Rev. A 74, 013607 (2006). 

V. M. Stojanovic, C. Wu, W. V. Liu, and S. Das Sarma, 

Phys. Rev. Lett. 101, 125301 (2008). 

C. Wu, Modern Physics Letters B 23, 1 (2009). 

A. Isacsson and S. M. Girvin, Phys. Rev. A 72, 053604 

(2005). 

A. B. Kuklov, Phys. Rev. Lett. 97, 110405 (2006). 

T. Miiller, S. Foiling, A. Widera, and I. Bloch, Phys. 

Rev. Lett. 99, 200405 (2007). 

G. Wirth, M. Olschlager, and A. Hemmerich, 
arXiv: 1006.0509 (2010). 

C. Wu, D. Bergman, L. Balents, and S. Das Sarma, Phys. 

Rev. Lett. 99, 070401 (2007). 

C. Wu, Phys. Rev. Lett. 100, 200406 (2008). 

C. Wu and S. Das Sarma, Phys. Rev. B 77, 235107 

(2008) . 

S. Zhang, H.-h. Hung, and C. Wu, Phys. Rev. A 82, 
053618 (2010). 

W.-C. Lee, C. Wu, and S. Das Sarma, Phys. Rev. A 82, 
053611 (2010). 

H. -h. Hung, W.-C. Lee, and C. Wu, arXiv:0910.0507 

(2009) . 

S.-L. Zhu, B. Wang, and L. M. Duan, Phys. Rev. Lett. 
98, 260402 (2007). 

K. L. Lee, B. Gremaud, R. Han, B.-G. Englert, and 
C. Miniatura, Phys. Rev. A 80, 043411 (2009). 
K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, ArXiv 
e-prints (2010), 1012.5864. 

T. Neupert, L. Santos, C. Chamon, and C. Mudry, ArXiv 
e-prints (2010), 1012.4723. 

E. Tang, J. Mei, and X. Wen, ArXiv e-prints (2010), 
1012.2930. 

A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
P. Bak, Rep. Prog. Phys. 45, 587 (1982). 
P. Bak, Phys. Today 39, 38 (1986). 

R. Jordens, N. Strohmaicr, K. Gunter, H. Moritz, and 

T. Esslinger, Nature 455, 204 (2008). 

J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, 

I. Bloch, and S. Kuhr, Nature 467, 68 (2010). 

U. Schneider, L. Hackermuller, S. Will, T. Best, I. Bloch, 
T. A. Costi, R. W. Helmes, D. Rasch, and A. Rosch, 
Science 322, 1520 (2008). 

W. S. Bakr, J. I. Gillen, A. Peng, S. Foelling, and 
M. Greiner, Nature 462, 74 (2009). 

R. O. Umucalilar, H. Zhai, and M. O. Oktel, Phys. rev. 
lett. 100, 070402 (2008). 

P. C. Haljan, I. Coddington, P. Engels, and E. A. Cornell, 



15 



Phys. Rev. lett. 87, 210403 (2001). [87] P. Streda, J. Phys. C 15, L717 (1982). 

[86] Y.-J. Lin, R. L. Compton, K. J. Garcia, J. V. Porto, and 
I. B. Spielman, Nature 462, 628 (2009). 



